Broken space-time symmetries and mechanisms of rectification of ac fields by 

nonlinear (non)adiabatic response 
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We consider low-dimensional dynamical systems exposed 
to a heat bath and to additional ac fields. The presence of 
these ac fields may lead to a breaking of certain spatial or 
temporal symmetries which in turn cause nonzero averages 
of relevant observables. Nonlinear (non)adiabatic response is 
employed to explain the effect. We consider a case of a particle 
in a periodic potential as an example and discuss the relevant 
symmetry breakings and the mechanisms of rectification of 
the current in such a system. 



I. INTRODUCTION 

Much has been written on noise induced transport, 
where noise may be coloured, or simply white plus 
time periodic signal. Arguments include violation of 
fluctuation-dissipation theorems, breaking of reflection 
symmetries of potentials in space, Maxwell daemons, 
mixing of harmonics (of e.g. a periodic drive), etc. All 
this applied to both classical and quantum systems, and 
extended at least conceptually to stochastic resonance, 
quantum stochastic resonance, etc. A recent review on 
ratchet transport by Reimann [Q provides with a lot of 
theoretical and experimental results, and we refer the in- 
terested reader to this work. 

For a set of related problems, such as directed particle 
current B0| , directed energy current M , average magne- 
tization |q]7 and nonlinear Hall effect jpj , to name a few, 
a recently published symmetry approach was shown to 
systematically account for all relevant symmetries which 
have to be broken in order to explain the observed recti- 
fication effect. The purpose of this paper is to generalize 
this approach, and to apply it to different physically rel- 
evant situations like underdamped, overdamped or zero 
temperature cases. We also use the nonlinear response 
arguments to make the results of symmetry considera- 
tions very transparent. We will argue for a rather minor 
role of additional fluctuations which are mainly responsi- 
ble for the area of phase space explored. Although we will 
discuss mostly the case of time-periodic external fields, 
we will also show how the nonlinear response concept can 
be generalized to fields which are quasipcriodic in time. 



The paper has the following structure. In the next 
section we discuss the symmetry breaking by consider- 
ing the adiabatic limit and making use of some general 



forms of nonlinear response functions. Section [II is de- 
voted to the detailed symmetry analysis of a classical 
particle moving in a periodic potential under the influ- 
ence of an external ac drive. Symmetry breaking leads 
to a directed current here. We will discuss the presence 
of ballistic channels which provide with a mechanism of 
current rectification. We will show that these ballistic 
channels survive in the presence of dissipation. Section 
IV addresses the case of a particle with nonparabolic dis- 
persion. The addition of a spatial potential allows us to 
tune the system in such a way that ballistic channels dis- 
appear. This will lead to an expected dramatic drop of 
the directed current value. Conclusions and discussions 
are given in Section ^ 



II. ADIABATIC RESPONSE TO SLOW 
AC FIELDS 

Consider a certain system in contact with a heat bath. 
The system is characterized by some internal (nonlinear) 
dynamics, and we will discuss concrete models below. So 
far we just need to know that its state can be charac- 
terized by certain variables which are functions of time. 
These could be functions of phase space variables of a 
classical system, or expectation values of operators for a 
quantum system. Consider one of such variables which 
we denote as A{t). The coupling to the heat bath will 
be characterized by at least two parameters - the tem- 
perature of the bath (we will use the notation of inverse 
temperature /?) and some set of relaxation times, in the 
simplest case just one relaxation time (we will use the no- 
tation of the inverse relaxation time i/). We also assume 
that the chosen variable A{t) is zero on average. 

Let us now apply a static field E which couples to A 
such that a nonzero average A = {A{t))t is generated. 
Its dependence on E is assumed to be given by some sin- 
gle valued response function A = f{E) (the physical 
meaning of the response function can be different, elec- 
tric polarization or magnetization, for example). The 
above mentioned symmetries will be connected with cor- 
responding symmetries of the function f{E). This func- 
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tion will be in general nonlinear, yet in some cases it may 
be expanded in a Taylor series around E = and start 
with a linear term (this term describes then the linear 
response). There are two possibilities: the function 
is either antisymmetric f{E) — —f{—E) or it is asym- 
metric /{E) + -fi-E). 

Before analyzing the case of adiabatically slow periodic 
ac fields, we will more rigorously introduce the notion 
of possible symmetries of an arbitrary periodic function 
with zero mean. 



A. Classification of symmetries of a periodic 
function with zero mean 

Consider a periodic function g(z + 2tt) — g{z) hav- 
ing zero mean. First such functions may be symmet- 
ric around certain z values. Without loss of generality 
this point may be set to zero and we find g{z) = g{—z). 
We will use the abbreviation in such a case. Second 
such functions may be antisymmetric g{z) = —g(—z) 
(abbrevation ga). Note that the points around which 
a function is symmetric and antisymmetric will be dif- 
ferent if the function possesses both symmetries. Fi- 
nally the function may posses shift symmetry (which is 
also called antiperiodicity) : g{z + t:) = —g{z) (abbre- 
viation gsh)- A zero mean periodic function may pos- 
sess none of the above symemtries, precisely one of them 
or all three simultaneously. In particular g^ functions 
may be expanded in a pure cosine Fourier series, ga 
functions in a pure sine Fourier series, and gsh func- 
tions show up with zero even Fourier components in 
their Fourier series expansion. As a consequence the 
simplest function g{z) = cos(z + zq) possesses all three 
symmetries. The function g{z) — cos(z) -I- cos(2z -I- zq) 
does not posses any of the listed symmetries except for 
Zq — 0,7r (gs) and zq = ±7r/2 (^q). A final example 
g{z) = cos(z) -I- cos(3z -I- zq) always possesses shift sym- 
metry (gsh) and in addition may be symmetric and anti- 
symmetric for Zq = 0,7r. 

Let us finally note that most cases under considera- 
tion use several harmonics contained in the drive. How- 
ever it may be also important to use pulse sequences. 
The symmetry considerations can be straightforwardly 
applied also to such cases. 

B. Periodic ac fields 

Let us now assume that the field E is slowly varying in 
time. Slowly means that the characteristic time scales of 
changes of E(t) are much larger than all other time scales 
in the system. Then we are dealing with the adiabatic 
limit of the response to time-dependent fields. We may 
use the response / to a static field and simply insert 
the slow time dependence: f[E{t)]. Consider a time- 
periodic field E{t) = E{t+T) with zero average E = 0. In 



general the corresponding value A = 1/T J f{E{t))dt will 
be nonzero for asymmetric /. However for antisymmetric 
single-valued / the average A will be zero provided the 
fimction E{t) is either antisymmetric E{t) = —E(—t) or 
has shift symmetry E{t) = -E{t + T/2). 
Consider as an example 

E{t) = El cos t + E2 cos(2t + a) + E3 sin(3i + a') . (1) 

For E2 — E3 — Eq. ([^) has all three mentioned sym- 
metries of periodic functions. For E^ = the function 
has no symmetries except for a = mr {n integer) where 
it is symmetric or for a = 7r(l/2 -I- n) where the function 
is antisymmetric. For E2 — Eq. (|^) has shift symmetry 
except for a' = 7r(l/2 -I- n) when it is in addition sym- 
metric and antisymmetric (around different time origins, 
of course) . 

Next we assume a case when the adiabatic response 
function / may be expanded in a Taylor series: 

f{E) = hE + hE'' + hE^ + ... (2) 

where the skipping of higher order terms is justified by 
the smallness of E. For /2 ^ the response is asymmet- 
ric and consequently any field from (|l]) will in general 
lead to a nonzero average of /. For /2 = however all 
shift symmetric and antisymmetric E{t) will yield zero 
averaged / (because it leads also to a shift symmetry or 
antisymmetry of / in time). As an example of mixing of 
harmonics the leading nonzero contribution for E2 ^ Q, 
i?3 = equals 

f [E{t)]dt^y3ElE2 cos a + ... . (3) 

In accordance with the above said / vanishes for a = 
7r(l/2-|-n). Note that the antisymmetry property of f{E) 
is linked to some internal symmetries of the dynamical 
system under consideration, as will be shown below. 

It is important to note that the written above is valid 
for single- valued response functions. If the coupling to 
a heat bath is too weak, response functions may contain 
hysteretic loops. Such cases need separate discussion 0. 

C. Quasiperiodic ac fields 

We may even consider quasiperiodic driving here. Sup- 
pose that we drive the system with a field E{t) — 
ei(t) -\- e2{t) where both components ei_2 are time pe- 
riodic, but with incommensurable periods. The resulting 
field E{t) would be then quasiperiodic. If considering low 
amplitude fields and expanding the response function to 
second or third order, the full dc average of A{t) will 
be given by a mere sum of the averages obtained in the 
presence of only one of the two periodic fields. Consider 
e.g. E{t) = ei cos Wit + 62 cos 0^2^ and an asymmetric re- 
sponse. In lowest order in the field amplitudes we obtain 
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(4) 

In higher orders in the response interference effects ap- 
pear, so that products of the field components enter the 
result. Yet for the lowest order case (/2) only the sum ap- 
pears. So we may assume that further changing the drive 
by adding more time periodic field components does not 
change the main result of a possibility of a nonzero aver- 
age. 

If we consider an antisymmetric response function /, 
we need to choose a more sophisticated quasiperiodic 
drive. Take e.g. E{t) = ei{t) + 62 (i) with ej(t) = 
en cosujit + ei2 cos(2w,;t -I- at). Assuming again that the 
field amplitude is small, we obtain in lowest order 

3 

/ = ^/a (eiiei2 cos ai -1-621622 cos ■ (5) 

Again we find that the periodic components of the 
quasiperiodic drive contribute additively in lowest order. 
Note that each contribution can be obtained by consider- 
ing a reduced drive function consisting only of this part. 
The corresponding contribution is clearly related to the 
symmetry properties of this reduced drive. 

D. Do we need more? 

Once a given problem is considered in the adiabatic 
limit like discussed above, the appearance of nonzero 
averages due to ac fields can be obtained, and further 
changing of parameters away from the adiabatic limit 
will change numbers, but not the fact of nonzero aver- 
ages. Ratchet transport in its general form falls into such 
a class of systems. Indeed a particle moving in a peri- 
odic potential V{x) = V{x + X) under the infiuence of an 
external colored noise and dissipation may be taken into 
the limit of zero temperature (no noise), slow periodic 
driving and zero mass (over damped case): 

x + V'{x)+E{t) = . (6) 

The response of (i) to a corresponding static field E is 
well known and has been discussed in connection with 
Josephson current-voltage characteristics (for symmetric 
potentials). The response is clearly nonlinear, and an- 
tisymmetric if the potential is symmetric in space. In 
case of asymmetric ratchet potentials the response is 
nonlinear but asymmetric. In the latter case a simple 
E{t) ^ cost signal with large enough amplitude will gen- 
erate a nonzero averaged current. In case of symmetric 
potentials a drive which is neither antisymmetric nor shift 
symmetric is needed. 

Nonadiabatic corrections are of interest and impor- 
tance. E.g. upon constantly increasing the frequency 
of a drive, currents (or other averages) may change sign. 



or increase or decrease by orders of magnitude [||, 
|lO| . Especially interesting are cases when certain (nona- 
diabatic) parameter limits give rise to new symmetries. 
These symmetries may be traced back away from the 
corresponding limits and help understand lots of pecu- 
liar features which are observed in such driven systems. 

If the adiabatic response simply vanishes, nonadiabatic 
response contributions will be not just mere corrections. 
In such cases the symmetry considerations of the under- 
lying equations of motion presented below appear to be 
the most direct way to predict rectification effects. One 
example of such a case is given in Q where the problem of 
driven quantum spins in ac magnetic fields is considered. 
While the adiabatic limit provides with zero induced 
magnetization component in y-direction there, nonadi- 
abatic response terms make this component nonzero, as 
expected from symmetry considerations. 

III. DIRECTED CURRENTS - THE CASE OF A 
CLASSICAL PARTICLE IN A PERIODIC 
POTENTIAL 

In order to make further progress in the understanding 
of rectification it is useful to define a model. Here we will 
consider the case of a particle of the unit mass moving 
in a space-periodic potential V{x) — V{x + A) under the 
influence of friction and external forces: 

mx + 'yx- f{x) - x{t) = , (7) 

where f{x) — —V'{x) and x{t) is zero on average. The 
mass m is assumed to be equal to one if not stated oth- 
erwise. If x{t) is a Gaussian white noise, the particle un- 
dertakes diffusion with zero net current, in accord with 
the fluctuation-dissipation theorem. If however x{t) con- 
tains correlations (color), then it is known that a nonzero 
net (dc) current is possible. In order to understand this 
result, it is appropriate to make the correlations in the 
noise as transparent as possible. The easiest way to do 
so is to choose 

X{t) - m + m , E{t) = E{t + T) , (8) 

where ^(t) is a Gaussian white noise and the time- 
periodic external field E{t) has zero mean. A proba- 
bilistic description of the system is then possible using a 
Fokker-Planck equation (see Appendix A) . 

In a next radical step we skip the f (t) term. The rea- 
son for that is that we will be left with a determinis- 
tic equation, whose symmetries may be studied. These 
symmetries involving operations in time t are assumed to 
hold even in the presence of a Gaussian white noise, since 
such a term does not contain temporal correlations. Note 
that the general nonintegrability of the resulting equa- 
tions may provide with irregular (chaotic) trajectories. 
Thus we are left with an equation of the form 

mx 7i - f{x) - E{t) = . (9) 
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A. The relevant symmetries and ways of violations 

In order to characterize the symmetries of Eq. (|^) we 
remind that the phase space dimension is three. As we 
look for nonzero average currents which are characterized 
by the velocity i, we have to check whether there exist 
symmetries which allow to generate out of a specific tra- 
jectory of (j^) another one with reversed velocities. The 
transformations of interest have to involve a change of 
sign of X leaving its absolute value unchanged. Thus we 
look for transformations which leave Eq. (m) invariant and 



either change the sign of x 
ously shift time t t + 



-X and simultane- 



• invert time t 



't and shift coordinate x 



X+Xq. 



The following symmetries can be identified 

T 

Sa : X ^ -x , t + — , if {fa , Esh} 



Si) '. X ^ X ^ t 



if {E, ,7 = 0} ; (10) 



Sc ■■ x^x+-,t^-t, if {fsh , Ea , m = 0} . 

Note that Sb and Sc require 7 = (Hamiltonian case) 
and m — (overdamped case) respectively. Another ob- 
servation is that all symmetries require certain symmetry 
properties of the external drive E{t), while the properties 
of the space-periodic force / does not matter for Sb- A 
proper choice of the drive may thus break all of the above 
listed symmetries with any coordinate dependence of the 
force f{x). 

It is useful to provide an interpretation of the action 
of the above symmetry operations on trajectories of (H). 
If the equation of motion is invariant under a symme- 
try, it implies that a given solution or trajectory, when 
transformed using the symmetry operation, yields again 
a solution or trajectory of the system. It may be the 
same trajectory or a different one. The symmetry Sa for 
7 7^ transforms an attractor into an attractor, repellers 
into repellers and basins of attraction into basins of at- 
traction. Thus if an attractor is mapped onto itself, the 
average velocity on that attractor is zero. If an attrac- 
tor is mapped onto another attractor, we find complete 
symmetry between the two attractors and their basins of 
attraction, while the average velocities will be of opposite 
sign. The symmetry Sb, which is valid for the Hamilto- 
nian case, simply relates two trajectories to each other, 
both having opposite average velocities. 

Once we add white noise, we have to consider some way 
of weighting different trajectories. In that sense, if two 
different trajectories are related by a symmetry Sa or Sb 
and have identical statistical weight, their contributions 
to a total average current will annihilate. 

The symmetry Sc is acting in the overdamped limit, 
yet it is reverting the sign of time. Thus it maps an at- 
tractor on a repeller and the basin of attraction of an at- 
tractor into the basin of repulsion of a repeller. It follows 



that due to the low dimension d = 2 of the phase space 
of this problem, the compactness of the phase space and 
the uniqueness of a trajectory running through a point in 
it, that all trajectories have zero average velocity when 
Sc is valid. 

In order to correctly incorporate the effects of noise, 
kinetic equations have to be considered. This is done in 
Appendix A. 

If at least one of the symmetries ( p^ ) is valid, we con- 
clude that the averaged current is zero if an additional 
white Gaussian noise term is added. This noise will sim- 
ply lead to an exploration of the whole phase space and 
thus realize on finite times different trajectories of the 
deterministic system. 

Once all of the above symmetries are broken, we 
may expect that in general the average current will be 
nonzero. The understanding of the mechanisms which 
will lead to this current can be best obtained in the ab- 
sence of noise for the deterministic system. 



B. Rectification mechanisms in the Hamiltonian 
limit 

As it was shown in |^ the value of the average current 
is strongly enhanced for the underdamped case, when 
the dissipation rate tends to zero. Thus we start our 
study of the mechanisms with the dissipationless case 
7 = 0. This Hamiltonian limit of Eq. (^) is generically 
characterized by nonintegrable dynamics with a mixed 
phase space containing both chaotic and regular areas 
. A stochastic layer appears around the separatrix of 
the integrable nondriven [E{t) = 0] Eq. (^. 

Contributions to a nonzero current may come from tra- 
jectories inside the stochastic layer as well as from regu- 
lar unbound motion. Because the latter always exist for 
motion in both directions, strong effective cancellation 
of currents takes place while the strongest remain- 
ing contribution again comes from the stochastic layer, 
which will be considered below. 

In the following we choose the following functions and 
parameters: 



V{x) — — cosx , 

E{t) = El cos ujt + E2 cos{2ujt + a) . 



(11) 



This choice ensures that the symmetries Sa and Sb are 
broken. In Figj^ we show the Poincare map of the main 
stochastic layer of our model. 
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FIG. 1. Poincare map of the main stochastic layer of (^) 
with 7 = and functions Q with Ei = 3.26, E2 = 1.2, 
a = n/2 and u = 1. The two filled circles correspond to the 
two observed limit cycle attractors in the weakly dissipative 
case with damping 7 — 10~*. 

Points with coordinates p = x and x mod 2t: are drawn 
after each period of the drive. In Fig. ^ the time depen- 
dence of the coordinate x(t) is shown. 
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in these boundaries and perform ballistic-like (regular- 
Uke) dynamics. The symmetry breaking of the equations 
of motion is reflected in a desymmetrization of these frac- 
tal boundary structures for the upper and lower bound- 
aries. This in turn leads to a desymmetrization of dis- 
tribution functions which characterize the probability to 
stick and stay in such a boundary region |l^] . In Fig.^(b) 
the probability distribution functions of the flight dura- 
tions to the left and the right (for technical details see 
Jl2t ) are shown. These functions are characterized by 
algebaically decaying tails, and we clearly observe the 
above mentioned desymmetrization. 
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FIG. 3. Probability distributions of ballistic flights to the 
right (solid line) and left (dashed line), (a) - case with dis- 
sipation and noise, (b) - Hamiltonian case. Inset: Length of 
ballistic flight versus time of ballistic flight for Hamiltonian 



In the inset of FigJ|(b) each point denotes the distance 
covered in a given flight during the time of ballistic-like 
motion. We observe a fine structure with three branches. 
Two major ones correspond to the main ballistic chan- 
nels with opposite velocities while and a minor third one 
corresponds to a channel with smaller negative velocity. 



C. Presence of dissipation 



FIG. 2. x{t) of trajectory from Fig.^. Inset: Zoom of x(t). 
Note the axis scales. 

We observe a drift in accordance with the symmetry 
breaking (|. The average velocity is approximately (i) « 
0.85. 

In order to understand the dynamical mechanisms of 
a nonzero dc current in the stochastic layer, we first note 
that the stochastic layer is bounded in phase space. The 
boundary contains of a fractal set of regular islands em- 
bedded in the stochastic layer. A trajectory from the 
stochastic layer may become trapped for quite long times 



When dissipation is present in the system (i.e. 7 7^ 0), 
the phase space of the system separates into basins of 
attraction of different low-dimensional attractors. Close 
to the Hamiltonian limit these attractors are limit cycles: 
x{t + T) = x{t) + 27rm, x{t + T) = x{t), me Z. These 
limit cycles are locked to the external periodic drive E{t), 
therefore their period T is characterized by T = 2mr/uj, 
n £ Z . The average velocity on the limit is given by 



, . , 1 / . m 

(x) = — / xdt = —uj 



(12) 



5 



When further away from the Hamiltonian Umit chaotic 
attractors can appear via period-doubUng bifurcations 

A numerical test reveals that in the case of 7 = 10 
only two limit cycle attractors appear. Their location in 
the Poincare map is shown in Fig. 0by two filled circles. 
The dependence of x(t) on both attractors is shown in 
Fig. @ 
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FIG. 4. p{t) = x{t) of both attractors for 7 = 10"* which 
correspond to the two filled circles in Fig. |l| 

Note that the attractors are located inside regular is- 
lands of the corresponding nondissipative system. These 
islands are characterized with nonzero winding numbers, 
and the sticking of the chaotic trajectory of the nondis- 
sipative case provides with the above discussed ballistic 
channels there. So we observe that the ballistic channels 
of the Hamiltonian system survive by transforming into 
limit cycle attractors of the weakly dissipative case. 

As the attractors are located inside the stochastic layer 
of the Hamiltonian limit, their basins of attraction are 
expected to show up with a complex folding structure. 
This fact is manifested beautifully if we add a noise of 
weak intensity which corresponds to a small temperature 
l//3 = 0.05: 



mm) = 2i/mt-t') 



(13) 



1//9 is small compared to the energy barrier of the pe- 
riodic potential. As it turns out, in this case of strong 
external driving the system is so far from the equilib- 
rium case, that different scales have to be used in or- 
der to compare with the noise intensity. In Fig. |^ we 
show the corresponding Poincare map over a total time 
of t/ = 10^. We observe that the trajectory is sticking 
for long times to the two attractors. But most impor- 
tantly we observe frequent escapes from the attractors 
basins. Once the trajectory is kicked out of such a basin, 
it starts to quickly explore the stochastic layer space, due 
to the weak damping, weak noise, and the above men- 



tioned complex folding of the basin boundaries. This 
results in the fact that the probability distribution of the 
velocities is far from being Maxwellian. 
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FIG. 5. Poincare map of the dissipative case with weak 
noise. The two filled circles correspond to the two observed 
limit cycle attractors for 7 = 10~*. 

The dependence of x{t) for the dissipative case with 
noise is shown in Fig. ^. Although the dc current value 
has changed compared to the Hamiltonian case, again the 
evolution is characterized by sticking to ballistic chan- 
nels. 
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FIG. 6. x{t) of trajectory from Fig.^ . Inset: Zoom of x{t). 
Note the axis scales. 

An evaluation of the corresponding ballistic flight time 
distributions is shown in Fig.|3|(a) . A remarkable similar- 
ity to the case of the distributions for the Hamiltonian 
case is observed. Also we observe that the power law tails 
of the Hamiltonian case are replaced by exponential ones 
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in the dissipative case with noise, due to an expected 
noise-induced cutoff in the maximum correlation time. 

A consequence of the above results is that the loss of 
ballistic channels due to the variation of some parame- 
ter may lead to a crossover-like decrease of the current. 
In the next section we will design a model which does 
possess these properties. 



IV. DIRECTED CURRENTS - THE CASE OF 
NONPARABOLIC DISPERSION 

So far we have discussed the case of a classical particle 
with parabolic dispersion, i.e. a kinetic energy quadratic 
in the momentum p. What happens if we consider non- 
parabolic dispersion? A prominent example would be a 
periodic dependence on p as for (— cosp), which reminds 
the consideration of a quantum particle evolution in one 
band approximations. 



A. Absence of additional potential 

Let us first consider the case of a particle with a gen- 
eral dispersion relation in the presence of additional ac 
driving: 



H = e{p)-xE{t) , e{p) = e{-p) 



(14) 



If the function e{p) is chosen to be periodic, its period is 
defined as Ap. Let us first consider the symmetries of the 
equations of motion which change the sign of the velocity 
x: 



i^e'ip) , p^Eit) . 
The following symmetries can be identified: 



(15) 



X - 


-X , p- 


-^-p,t^t+- if {Esh} , 


(16) 


t - 


' -t , p- 


> if {Es} , 


(17) 


X - 


-X , p- 


+ Y if Wsh} > 


(18) 


t - 


' -t , p- 


>P+^ if U'sH,K}- 


(19) 



Note that the last two operations (jl^JT^) may apply only 
for periodic e(p) functions. Furthermore these operations 
change the energy of the undrivcn system. If these op- 
erations connect different trajectories, they should not 
matter at finite temperatures since different energies con- 
tribute with different weights. 

Choosing the free particle case e(p) = we may 

arrive at the conclusion that both relevant symmetries 
( |l6| , p7|) can be violated by a proper choice of E{t). Yet, 
as shown in ||^,|l^, the expected dc current should be 
zero in such a case, in fact independent on the strength 
of some additionally applied dissipation and noise. The 



case of e(p) = — cosp is more involved. On one side it is 
well known that the kinetic Boltzmann equation solution 
provides with a nonzero dc current if both symmetries 
( |l6| , ^7|) are broken On the other side this dc current 
tends to zero as the dissipationlcss Hamiltonian case is 
approached. We will show in the following that the rea- 
son for that is an additional symmetry of the solutions 
of (|l|) due to the integrability of (|l^ . 
Indeed, integration of ( p^ ) yields 

p{t) = g{t;po,to) = PQ + Ei{t) - Ei{to) . (20) 

where Ei{t) — E{t). The function g has the property 

9{x;y,z) ''g{z;-y,x) . (21) 

This symmetry operation relates points on a given tra- 
jectory to points on a set of other trajectories. Such a 
symmetry can not be derived from the equations of mo- 
tion. It rather is a result of the integrability and reflects 
the symmetry of the solutions, i.e. of the phase space flow 
(in contrast symmetries of equations of motion relate two 
trajectories to each other). 

The symmetry (21) is in place independent of the 
choice of the functions e{p) and E{t) and provides with 
the following consequences. If we consider a correspond- 
ing kinetic equation with finite dissipation, the loss of 
correlations implies an averaging over the initial phase 
to of the field E{t). This averaging persists in the dissi- 
pationlcss limit, while on the other side we recover the 
above considered Hamiltonian properties. The averaging 
over to leads to an enforcement of the symmetry (|2l|), 
which in turn by changing the sign of p changes the sign 
of X. Thus a vanishing dc current is expected in that 
limit, in accordance with numerical and analytical analy- 
sis of kinetic equations. This result can be also obtained 
from (|0l) and (|l|) if an averaging with a distribution 
function, which is symmetric in po and independent of x 
and is performed: 



dt 



dtox 



po 



dt 



dtnxl 



-po 



0. 



(22) 



To detect a dc current carried by the electrons in a sin- 
gle band one has to break the integrability of the semi- 
classical equations of motion. We consider this case in the 
next section. An alternative way is a much more careful 
treatment of the quantum mechanical problem (see, for 
example |[T5[) which is however far beyond the scope of 
the present paper. 



B. Presence of additional potential 

Things change drastically if we add an external space- 
periodic potential to (p^). Such a system is in general 
nonintegrable, and we may expect the additional sym- 
metry ( [2l| ) to not be present. To be precise, we choose 
the following model here: 
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1 ^2 ( X 
H ^ — cosp cos(2iD) cos — 

2 ^ 3 



xE{t) . (23) 



The corresponding equations of motion read 



sinp + sin(2p) 



(24) 
(25) 



The free parameter ^ can be used to reach the hmit of 
(^. Indeed for ^ ^ p 0. Choosing a small value 
of \p\ <^ 1 will keep the momentum small. This allows 
for an expansion of the right hand side of (|^ to first 
order. Additional rescaling x ^ x/£, and p {'^/Op 
transforms the problem exactly to Eq. (||) (with 7 = 0). 
On the other side, a choice of a finite value of p may 
lead to trajectories which are not contained in Eq. (H). 
In addition (23) is periodic in p, so the phase space is 
compact. 

Let us discuss the relevant symmetries of Eqs. ( P4[) - 
(p5|). Due to the presence of the cos(2p) term in ( |23| ) 
the shift symmetry of the kinetic energy is broken (m 
terms of electrons this implies loosing of the particle- 
hole symmetry). We can identify the following symmetry 
operations which lead to a change of sign of x: 



X —X , p 
t-> -t , p- 



-p , t 
-P if {E,} 



if {E,h}; (26) 



(27) 



In the following we choose E(t) from ( pT[ ) which ensures 
that both symmetries are broken. 

For C ^ 1 the phase space of ( p3| ) is characterized by 
the presence of two disconnected stochastic layers. One 
of them corresponds to the layer in Fig. |l|. 



continuously transforms into the layer in Fig. |l| when 
^ ^ 0. An increasing of ^ will lead to a merging of 
the layers, which is followed by losses (closings) of ballis- 
tic channels, that are located near the boundaries of the 
central stochastic layer. The Poincare map for £, = I 
shown in Fig. S. 



is 



D. 



-4 







FIG. 8. Poincare map of Eqs. (§^-(||) for f = 1. 

Let us discuss a trajectory with initial conditions p = 
and X — ^TT, which ensures that we always start in the 
stochastic layer which corresponds to Fig.l in the limit 
of small ^. The dependence of the scaled dc current js — 
is shown in Fig. || as a function of ^. Specifically 
its values for the presented Poincare maps are js ~ 0.63 
for ^ = 0.35 and js ~ 0.002 for ^ = 1. We observe indeed 
that the current variations correspond to a closing or 
opening of ballistic channels. 



D- 







FIG. 7. Poincare map of Eqs. (M)-(N) for ^ = 0.35. 



In Fig. ^ we show the Poincare map of (23) for 
^ — 0.35. Note that the central layer is the one which 
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FIG. 9. js as a function of ^. Note that the error in de- 
termining the current value is of the order of 10%. A fine 
structure of pronounced peaks is connected with the opening 
of new ballistic channels due to the overlap of the stochastic 
layer with high-order resonances. 

The above results support the general expectation that 
in the Hamiltonian hmit a mixed phase space is needed, 
so that the current vanishes both in the case of an in- 
tegrable system as well as in the case of a fully chaotic 
one. 



V. CONCLUSIONS 

We presented a symmetry approach to the effect of 
rectification due to external ac fields applied to a low- 
dimensional dynamical system with optional contact to 
a heat bath. The nonlinear (nonadiabatic) response is 
used to explain the effect. We explained the mechanisms 
of such a rectification for different cases. In the case of 
a particle moving in a space-periodic potential the ex- 
planation is given in terms of desymmetrization of bal- 
listic channels which correspond to motion in different 
directions. While this explanation starts from the case 
of a nondissipative Hamiltonian limit, we showed that 
the ballistic channels are robust with respect to adding 
dissipation and fluctuations. The mechanism of the di- 
rected current in the presence of dissipation is hidden 
in the desymmetrization of the basins of attraction of 
previously symmetry-related limit cycles with oppositely 
directed velocities. 

A recent geometric approach by H. Schanz et al in fl^ j 
provides with an elegant way to account for the average 
drift velocity in the Hamiltonian case. The basic ingre- 
dient of this approach is the assumption of an invariant 
density distribution in a stochastic layer being constant 
inside the layer. The resulting sum rule, when evaluated 
for concrete cases, provides with average velocities which 
quantitatively agree with our observations. A conse- 
quence of the sum rule is that a nonzero directed current 
may appear only in systems with a mixed phase space. 
This result correlates with our discussion of the mecha- 
nisms of directed currents in terms of ballistic channels, 
which may appear only provided the stochastic layer has 
some bounds with regular phase space regions, and the 
ballistic channels occur precisely in the neighborhood of 
these bounds. Open questions are: i) is there purely 
diffusive (chaotic) directed motion besides of transport 
through ballistic channels? (ii) when can one of these 
two mechanisms prevail? 

It was shown that the symmetry properties of the 
Fokker-Planck equation match the symmetries of the cor- 
responding deterministic equations of motion in the ab- 
sence of noise. Similar results hold for a Boltzmann equa- 
tion §. 

Our symmetry considerations may be as well used to 
explain rectifications in such diverse situations as i) the 



directed motion of particles in non-Newtonian liquids 
|l7| , ii) the appearance of ring currents for particles mov- 
ing in a two-dimensional space-periodic potential ||l8|,|l9| , 
iii) the appearance of directed heat fiux currents in sys- 
tems of interacting particles iv) the appearance of 
a nonzero dc polarization, v) quantum ratchets to 
name a few. Symmetry analysis is also instructive in the 
case of less conventional models, like systems with built- 
in asymmetry p0|-|2^ , where the asymmetry is hidden in 
a many body system's internal interactions. 
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Appendix A 

The Fokker-Planck equation for the probability distri- 
bution W{x,p,t) with p = i of equations (0),(H) reads 

i 



dW 
~dt 



dW 
dx 



m dp 



7 d 
mf3 dp 



W 



(28) 



This equation is linear in W, preserves the norm 
/ Wdx dp and is dissipative. For a fixed norm any initial 
condition will converge to a single attractor solution Wg ■ 
For the case E{t) = it is easy to see that the attractor 
is the Gibbs distribution. For nonzero E{t) the attrac- 
tor solution Ws will be periodic in x and t. The average 
current is given by 



(P) 



pWgdx dp dt 



(29) 



It follows that ( Pq ) is invariant under the following trans- 
formations: 



Sa : X 
Sb ■■ X ' 



■ -X , p 
X , t 



-P, t^i+l , if {/a , E,h} 



if , 7 = 0} . (30) 



At the same time ( [2^ ) is changing sign. Since the solution 
Ws is unique, the conclusion is that (p) = 0. 

The overdamped limit m = is described by the fol- 
lowing Fokker-Planck equation 

dW d 1 d'^W 

The average current is given by 

(p) = / f{x)Wsdx dt . (32) 
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The symmetry Sa is again holding. However the symme- 
try 



X + 



A 



t 



-t 



if {fsh , Ea , m ^ 0} 



does not follow from (|31 



in a similar way. It is interesting 



to note that Reimann [23[ considered the symmetries of 
the original stochastic differential equation and 
argued for a holding of Sc for nonzero noise intensities. 
Fistul showed that Sc connects the original Fokker- 
Planck equation ( |3l] ) with its conjugated counterpart. 

Here we will use a related approach and prove the van- 
ishing of (^) when Sc holds. First we remind that the 
distribution Ws is periodic in x and t. Next we note that 
the operators d/dt and d/dx are anti-Hermitian on the 
space of a;, t-periodic functions. Let us define the opera- 
tor T as 



T = 



dx 



mi, 



(33) 



Then (31) can be rewritten as a Lippmann-Schwinger 
type integral equation: 



W = l + Tf{x)W 



(34) 



Provided the conditions for Sc hold, i.e. E{t) = Ea{t), 
the operator T has the property 

T^f = -T{-t) , T{x + xo) --T{x) . (35) 

Expanding (|34|) in a formal series in f{x) we obtain 



0.04 



0.02 
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FIG. 10. (p) as a function of a (see text). 

For the above choice of functions f{x) and E{t) 
the symmetry Sa is violated because E{t) is not shift- 
symmetric. At the same time the symmetry Sc is violated 
for all values of a except for a = 0, ±7r, ±2tt, .... In Fig. 

the dependence of the average current {p) (^2|) on a 
is shown. We observe that for a = 0, tt, 2t: the symmetry 
Sc is restored and the current is vanishing (the absolute 
numerical value is less than 10^^^). 



71=1 



With (B2p we obtain for the average current 



(p) = T.j fi^){Tf{x)rdxdt 



(36) 



(37) 



Since all terms in (|3^) are real valued, we conclude that 
all integrals with even n vanish because of f{x) = fsh{x) 
and these terms contain odd powers of /. All integrals 
with odd n vanish because of (^), since they contain 
odd powers of T. Thus we find that indeed the average 
current exactly vanishes when Sc holds. 

We also performed a numerical solution of (^) for 
T = 7 = 1 and f{x) = sin.T, E{t) = sin(tjt)+sin(2wt-|-a) 
with Lo — 0.8. The function W was expanded in a Fourier 
series in x and twenty Fourier harmonics have been taken 
into account. The resulting coupled ordinary differen- 
tial equations were simply integrated until the system 
reached the final attractor W... 
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